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We study scalar mode perturbations (magnetosonic waves) induced by a helical stochastic cos- 
mological magnetic field and derive analytically the corresponding cosmic microwave background 
(CMB) temperature and polarization anisotropy angular power spectra. We show that the presence 
of a stochastic magnetic field, or an homogeneous magnetic field, influences the acoustic oscillation 
pattern of the CMB anisotropy power spectrum, effectively acting as a reduction of the baryon 
fraction. We find that the scalar magnetic energy density perturbation contribution to the CMB 
temperature anisotropy is small compared to the contribution to the CMB iS-polarization anisotropy. 

PACS numbers: 98.70.Vc, 98.80.-k 
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O ! I. INTRODUCTION 

(N , 

^ . A promising explanation for observed uniform magnetic fields in galaxies is that they are the amplified remnants of 

Qa seed cosmological magnetic field (for reviews see Refs. Jj]) generated in the early Universe .2, 3]. A seed magnetic 
field may have an helical part Magnetic helicity plays an important role in magnetohydrodynamical processes in 
■ the primordial plasma as well as in cosmological perturbation dynamics; In particular, magnetic helicity influences 
\ the inverse cascade mechanism — when energy is transferred from small to large scales — [5| , and as a result affects 
large-scale magnetic field formation Q. 

The average energy density and helicity of the magnetic field must be small to be consistent with the observed large- 
scale spatial isotropy of the Universe. In this case the linear theory of gravitational instability can be used to stud 
perturbation dynamics 0)11, B [13] ■ ^ cosmological magnetic field induces scalar, vector and tensor perturbations 
, At linear order each mode evolves separately. (At second order the modes are coupled and this results in non-Gaussian 
^ ■ effects [ni.) 

The vector (vorticity) and tensor ( grav itational waves) perturbation modes induced by a cosmological magnetic 
^ field have attracted a lot of interest [H, El [H, El [H, El, El, llSl- This is partially because they give rise 
to a i3-polarization CMB anisotropy signal, which vanishes for density (scalar mode) perturbations at linear order. 
f~i ■ Any cosmological signature of a primordial magnetic field is a potential probe (for a short review see Ref. [2l|). 
Qh' For example, the limit on a chemical-potential-like distortion of the CMB Planck spectrum leads to a limit on 
Q , a cosmological magnetic field of order 10~^ — 10~^ Gauss on 1 — 500 kpc length-scales Similar limits on a 

$H ' cosmological magnetic field generated during inflation Q are obtained from CMB temperature and polarization 
c/3 . anisotropy and non-gaussianity data from vector and tensor perturbation modes jl^ . 

Ci On the other hand, the scalar mode of magnetically driven perturbations also has a significant effect on CMB 

^ , fiuctuations: fast magnetosonic waves shift the CMB temperature anisotropy power spectrum acoustic peaks [H, [25l . 
12611 . In this paper we present a si/Btematic treatment of scalar magnetized perturbations that complements earlier work 
[Tol . El, m, m, m, [13, [28|. [ 291 [30l | . Using the total angular momentum formalism fsi'l and analytical approximation 
^ . techniques described in jl4| we obtain analytical expressions for CMB temperature and polarization anisotropies. 
■ ■ ' Reference [l^l presents numerical computations of scalar CMB temperature and polarization anisotropies in the 
case of a stochastic cosmological magnetic field with magnetic field power spectral indices ns = 1,2,3. Here we 
consider a general cosmological magnetic field and contrary to Ref. [27] we account for the Lorentz force term in the 
Euler equation for baryons, in accord with the analyses of Refs. [oi llOl. Il9l|. The main new results are approximate 
analytical expressions for the CMB fluctuations. This analysis allows us to identify two different effects arising from 
the magnetic field: (i) a rescaling of the photon-baryon fluid sound speed (that is responsible for the shift of the CMB 
acoustic peaks); and, (ii) effects from non-zero magnetic anisotropic stress (that is responsible for the additional CMB 
-E-polarization anisotropy signal). 
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The rest of this paper is organized as follows. In the next Section we describe the magnetic field, including the power 
spectrum, the anisotropic stress tensor, and its connection with the scalar (longitudial) part of the Lorentz force. We 
present analytical expressions for two-point correlations functions of the magnetic field energy density, the scalar 
part of the Lorentz force, and the magnetic field anisotropic stress. In Sec. Ill we derive the equations that govern 
scalar magnetic perturbations (the Einstein and matter conservation equations) and discuss general solutions and the 
initial conditions we adopt. In Sees. IV — VI we use semi-analytical approximations to compute CMB temperature 
and polarization anisotropics (as well as cross-correlations between temperature and ^^-polarization anisotropics). We 
conclude in Sec. VII. 



II. MAGNETIC FIELD STATISTICAL PROPERTIES 



A. Power spectrum 



We assume the presence of a Gaussianly-distributed stochastic helical cosmological magnetic field generated during 
or prior to the radiation-dominated epoch, with the energy density of the field a first-order perturbation to the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) homogeneous cosmological spacetime model. We neglect fiuid back- 
reaction onto the magnetic field, therefore the spatial and temporal dependence of the field separates, B(t,x) = 
B(x)/a^. Here a{t) is the cosmological scale factor, normalized to unity at the present time and B(x) is the 
magnetic field at the present time. Since the magnetic field energy density is first order, the magnetic field is 1/2 
order. 

Smoothing on a comoving length A with a Gaussian smoothing kernel oc exp[— x^/A^], we obtain the smoothed 
magnetic field with mean squared magnetic field — (B(x) • B(x))|a and mean squared magnetic helicity — 
A|(B(x) • [V X B(x)])|a jUE^I- Corresponding to the smoothing length A is the smoothing wavenumber k\ — 2n/\. 
We use 

B,ik) = I d\e^ -B,{^), B,(x) = j ^e-'^-i3,(k), (1) 

when Fourier transforming between position and wavenumber spaces. We assume flat spatial hypersurfaces (consistent 
with current observational indications, e.g., Ref. [33|). 

We also assume that the primordial plasma is a perfect conductor on all scales larger than the Silk damping 
wavelength As (the thickness of the last scattering surface) set by photon and neutrino diffusion. On much smaller 
scales we model magnetic field dampin g by an ultraviolet cut-off wavenumber = 271 /Xjj that is due to the damping 
of Alfven waves from photon viscosity [23, HOl (see Eq. (1) of Ref. [Ts^l); here Ad ^ Ag. 

Under these assumptions the magnetic field two-point correlation function in wavenumber space is 

(i?r^(k)B,(k')) = (27r)35(3)(k - k')[/^,,(k)PB(fc) + ie.,,ikPH{k)]. (2) 

Here i and j are spatial indices, i,j € (1,2,3), ki — ki/k a unit wavevector, Py (k) = 5ij — kikj the transverse 
plane projector, e,,,; the antisymmetric symbol, and (5''^^(k — k') the Dirac delta function. The power spectra of the 
symmetric and helical parts of the magnetic field, Psik) and Puik), are assumed to be simple power laws on large 
scales, 

2_-2\3r>2 n 2\:ijT2 

P5(fc)=P5ofc"- =^^^^(Afc)"-, p^[k) = PH,k-- = -——^(\kr- , k<kD, (3) 

and vanish on small scales where k > k]j. The spectral indexes Ub and Uh are constrained by the requirement of 
finiteness of mean magnetic field energy density {ns > —3) and mean magnetic helicity {nn > —4) in the infrared 
region at small k. In addition, causality requires Psik) > |Pff(A:)| (the Schwartz inequality), [33l |. 



B. Anisotropic stress-energy tensor 



We consider the effects of a cosmological magnetic field on the CMB at high redshift when the Universe is hot and 
the plasma a good conductor. As a result the magnetic field lines arc dragged by the matter fiuid and this generates 
a weak "frozen" electric field E — — v x B, where v is the perturbed (first order) 3-velocity of the fluid. We neglect 
this weak electric field in what follows since the energy density of this electric field contributes at third order in the 



perturbation expansion. At the current time, the space-space part of MaxweU stress-energy tensor for the magnetic 
field is 
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B.(x)i?,(x)-i5,,i?2(x) 



(4) 



where rjo is the current value of conformal time 77 — dt/a[t). The energy density and the anisotropic trace- free part 
of the space-space components of the stress-energy tensor of the magnetic field at the current time are 



Pb(x,?7o) 
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i32(x)--r''-(x,%), 



^/f (x,%) - Is^jT^' {x,f]o), 



(5) 
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where r'"' — 5ijT-^\ Both pb(x,77o) = pB{'x.,'n)a^ and t^^j^^(x, 770) — T-^\x,r])a'^ are quadratic in the magnetic field, 
and their wavenumber space transforms are convolutions of the magnetic field, 

d3p 



PS (k, 770) 



-r-(k,,„) = i-| ^i?,(p)B,(k-p), 



B,{p)B,{k - p) - ^S,jBi{p)B,{k - p) 



(7) 
(8) 



We assume that Eq. ([6]) is a first order perturbation and decompose it into scalar, vector, and tensor parts, 

T^^\x.) — H^pix.) + u!^\-x.) + n^J''(x) 0, HI]. This decomposition is more conveniently done in wavenumber space. 

We use projection operators and find n-J^(k) = [_Pib(k)fcj -I- Pjb(k.)ki]kaTabi^) (for the vector part) [13, [ll| and 

n,(J'(k, Tjo) = [P^a(k)Pj6(k) - (k)Fa6(k)/2]Tafc(k, Tjo) (for the tensor part) The scalar parts of pb and r^f ^ 

determine the scalar part of the magnetic source. The scalar part n|J^(k) has to be proportional to kikj — %/3 [3], 
so we define 



n!f(k)^^fM,-k,V^"Hk) 



(9) 



Here the scalar n("^''(k, 77) — n('^)(k, ?7o)/a'* is associated with the anisotropic stress of the magnetic field. (At leading 
order the isotropic pressure has the same time dependence, pb(x, 77) — pb(x, 770)70"*, and is related to the magnetic 
field energy density by pb(x, r;) = ps(x, 77)/3.) It is straightforward to determine n'^'^)(k) by applying knkm — 5nm/?> 
on rlm^(k), i.e.. 



n^'^Hk) — knkrnTnmO^J 



(10) 



where we use SnmTnm — 0. Our n''^)(k, 770) is related to the T''^)(k, 779) of Ref. ^ through r('^^(k, t^q) = 
3n(^)(k,77o)/2. 

The scalar n^'^^ is related to the scalar part of the Lorentz force L(x, t^q) = — [B(x) x (V x B(x))]/(47r) and the 

isotropic pressure ps(x, 770). We introduce a scalar L^^^ defined by l|'^'(x, 770) — iL^^\-x.,riQ), where hf^ is the 
scalar part of the Lorentz force. Using the Maxwell equation V • B = 0, the Lorentz force is 



i,(x,77o) = — 



S,(x)V,i?.(x)--V,i32(x) 



and the corresponding scalar part is derived through the scalar 



V2L(^^)(x,7;o) = V.Lf)(x,77o) = ^ 



1. 



(V.i3,(x))V,i3.(x)--V^i?^(x) 



where = ViV; is the Laplace operator and we have used V • B = 0. In position space Eq. PT!]) reads 

1. 



v2n(^)(x, 
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V,V,(B,(x)i?,(x))--V^B^(x 



(11) 



(12) 



(13) 



and comparison with Eq. p2p results in 

n(^)(x,,o) = ^^^+i(^)(x,^o), (14) 

the analog of Eq. (5.44) of Ref. [iG]. 

Since we consider a stochastic magnetic field we present here various correlations and averages of the magnetic 
source for scalar perturbations. The wavcnumber-space scalar two-point correlation function is 

(n(^)*(k,%)n(^)(k',,7o)) = {27rf\n^^\k,7jo)\'s^'\k-k'), (15) 

where the power spectrum |n^'^^(fc, ?7n)P depends only on fc = |k|. (The vector and tensor two-point correlation 
functions are given in Refs. [iJ, H^, [l8'| . ) Two-point correlation functions of the magnetic field energy density and 
the scalar part of the Lorentz force are defined in a similar manner, 

{p*s{Kijo)PB{k',vo)) = {2TTf\pB{k,m)\'S^'H^-^'), (16) 
(L(^)*(k,,7o)i(^Hk',r?o)) = (2^)3|L(S)(fc,^„)|25(3)(k_k'). (I7) 

Note that L^^\-K,r]) is related to J^(x, 77) of Ref. [i| through V2l('5)(x, 77) = J^(x, 7])/(47r), so L^'^^k, tj) = 
-^(k,,7)/(4^). 

The scalar power spectra |n'^'^)(A;)p, \pB{k)\^, and \L^^'>{k)\'^ are determined by the symmetric part of the magnetic 
field power spectrum Psik) and do not depend on the magnetic hehcity spectrum Pnik) (also see Ref. fll|). 

in^^Hfc)!' - ^ 1 d^p FB(p)Pi3(|k -p|) [9(1 -7')(1 -/?')- 6(1 7' -/3') + (1 + m')], (18) 

|PB(fc)P = 1^ J d^PBip)PB{\k-p\)il + p^), (19) 

\L^'\k)\' - ^ / d'p PB{p)PB{\k - p|) [47/5(7/? - M) + (1 + m')] • (20) 

Here 7 = k-p, /? = k-(k — p)/(fc|k — p|), and /i = p • (k — p)/(p|k — p|). The relations to the two-point correlation 
functions given in Ref. [U.] are (V*(k)T(k')) = {p*g{k)pB{k')) and (T*^(k)TS(k')) = 9(n(^)*(k)n(^) (k'))/4; T(k) and 
r^(k) are given in Eqs. (2.15) of Ref. [11[. 

Using the power law magnetic field power spectrum of Eq. Q, we can obtain expressions for the power spectra 
|n'^'^)(fc)P, |ps(fc)P, and |L('^'(fc)p in the semi-analytical approximation where we divide the integration range into 
p < k and p > k parts and consider two limiting ranges p <^ k and p ^ k [3, [l^ [l^. In particular, the magnetic 
energy density correlation power spectrum is^ 



\PB{k,rio)\ - 



32(2^5 + 3)r2(nB/2 + 3/2) 
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Ub I k 



riB + 3 \k 



(25) 



^ The power spectrum \pB(k)\^ is related to the two-point correlation function of the energy density of the magnetic field. By definition 
the r.m.s. magnetic field energy density ps''™'^ is 



1 



[ps™=(w)]' = (pb(x)ps(x)> = — / dkk'\pB{k,rjor- (21) 



2n' 

Replacing the upper limit of integration by the cut-off scale ku and using Eq. ^ below, we find 

PR™=(r)n) = — ' . (22) 

8V2(n-h3)7rr(nfl/2-h3/2) 

Pb™'^{vo) differs from the average magnetic energy density j5g = (ps) which is determined by the power spectrum EB{k) = k-^PB{k) /tt^ , 

1 1 r 

PbM = —{Bi{^)Bi{^)) = — dkEeik), (23) 



which, using the cut-off scale ko, gives 

B2(A:dA)"b+3 
An{nB + 3)r(ns/2 -|- 3/2) 

for ub = —3, ps(j7o) = B^/(8tt). We note, in particular, that the average Lorentz force vanishes but the r.m.s. Lorentz force is not 



For riB > — 3/2 this expression is dominated by the cut-ofF scale ku, and for large kjj it does not depend on k, while 
for UB < 3/2 we get \pBik,r]o)? oc fc2"B+3^ j^- ^e shown [H that |pB(fc,?7o)P - 9|L(^Hfc, %)l VS ~ 
and |pB(fc,ryo)P~9|n(^)(fc,r;o)|V4.' 

Increasing magnetic helicity reduces the vector part of Lorentz force two-point correlation function (L^^^'l^^'') ([^, 
[isj : this reduces parity-even CMB fluctuations E^E^); but leaves the scalar part unchanged (iL^'^)p is independent 
of magnetic helicity). In contrast, Refs. [ol. fl0l.l27j neglect the Lorentz force for a maximally helical magnetic field. 
They argue as follows: since (B • (V x B)) is maximal, the average of the Lorentz force oc (B x (V x B)) for a such 
field is minimal or even zero — this is a valid approximation for a homogeneous field, and results in the force-free 
approximation [§] — but this is not applicable for a stochastic field. In the case of a stochastic field the average 
Lorentz force is zero (as is the average magnetic field itself) but the Lorentz force two-point correlation is non-zero 
(see footnote 2 above). This affects stochastic peculiar motions (vorticity perturbations) of charged particles TS], and, 
as we show below, the dynamics of density perturbations also. In the stochastic field case the force-free approximation 
should be used with caution. 



III. SCALAR MAGNETIC PERTURBATIONS 



In this section we study the dynamics of linear magnetic energy density perturbations about a spatially-fiat FLRW 
background with scalar metric fluctuations. The metric tensor can be decomposed into a spatially homogeneous 
background part and a perturbation part, g^^ — gj^J + Sgf^,^; Greek letters are used for spacetime indices, fi,^ € 

(0, 1, 2, 3). For a spatially-flat model, and working with conformal time, the background FLRW metric g^j^J = a^rj^^, 
where -q^i, — diag(— 1, 1, 1, 1) is the Minkowski metric tensor. Scalar perturbations are gauge dependent because the 
mapping of coordinates between the perturbed physical manifold and the background is not unique. We work in 
Newtonian (longitudinal) gauge in which the metric tensor is shear free . Scalar perturbations to the geometry are 
then described by two scalar gravitational potentials 5* and $ where 

6gil^ = -2aH, Sg^p = "la^Mi,. (26) 

In this gauge, neglecting vector and tensor fluctuations, the line element ds^ = a^(r?)[— (l-|-25')d77^+(5ij(l + 2$)da;'da;^]. 
Matter perturbations are described by the perturbation of the complete stress-energy tensor, (5t^, which is the 

sum of the perturbed fluid and electromagnetic stress-energy tensors, br^ — '^jStJ^ + t^^^'^ . The subscript / 
denotes the three different matter components we consider, photons (7), baryons (6), or cold dark matter (c), which 
we model as fluids. A magnetic field infiuences perturbations without changing the background metric which is 
determined by the photon, baryon, and cold dark matter background densities. For simplicity we ignore neutrinos; 
Refs. 0, [13, [2^ account for the effects of these relativistic weakly interacting particles. Since we focus on dynamics 
at large redshift we also ignore a possible cosmological constant or dark energy. We decompose perturbations into 
plane waves oc exp(i k • x) and in what follows equations are presented in wavenumber space. 

The magnetic field source affects the motion of baryons. Before recombination the photon-baryon plasma is dom- 
inated by photons and has the relativistic equation of state p — pji. The baryon-photon momentum density ratio 
R(r\) = 3pb{'r])/{4p^{r])) — SaflQi, / (AftQ^) , where p-y and pb are the energy densities of photons and baryons, flo-y and 
flob are the photon and baryon density parameters measured today, i.e., Slg/ — pf{'r]o)/pci, where pcr = SttG / {3Hq) 
is the critical Einstein-de Sitter density and Ho is the Hubble constant. At early times i? <C 1; at the last scattering 
surface when photons and baryons decouple i?dec — 0.35. 

We consider a primordial magnetic field and assume that the stress-energy of the magnetic field is not compensated 
by anisotropic stress in the fluid. That is, we have non-zero initial gravitational potentials but vanishing initial fluid 
energy density and velocity perturbations. This is possible for a magnetic field generated during inflation [3, HI]. This 
will limit the spectral index of the magnetic field, ub < 0.'^ 

In subsections A and B below we present the linear perturbation theory equations for the metric and matter 
perturbations. In subsection C we discuss the speed of sound rescaling in the presence of a magnetic field, while in 
subsection D we consider initial conditions. 



2 It may be shown that L(^) ~ -272^^/3 ~ -ps and n(^) ~ -2pfl/3. 

^ This hmit on the magnetic field spectral index was obtained for a magnetic field generated through coupling of the inflation and 
hypercharge during inflation Ref. [sSl argues for the same result while Ref. [ssll relaxes the limit to rig < 2. 



A. Metric perturbations 



In zero-shear (Newtonian) gauge, part of Einstein's equation for scalar metric perturbations become two Poisson 
equations 0, @| , which in wavenumber space are 



PB 



(27) 



(28) 



Here (5/ is the /-th fluid density perturbation, is the scalar part of the /-th fluid velocity perturbation ' — kw^ , 

IIj is the scalar part of the anisotropic /-th fluid stress-energy tensor, 11^ — {kikj — Sij /3)Tf,ij , and an overdot 
represents a conformal time derivative d/drj. For economy of notation we do not explicitly show the wavevector (k) 
or time {rf} dependence of the variables in these and following equations. (For instance, the magnetic field energy 
density psiv) = Pb{i1o) I o-'^ ^ the Lorentz force L^^\r]) = L^^\rio)/a^, and the anisotropic stress of the magnetic field 
Ui^^ri) = n('5)(7?o)/a''-) Note that the combination on the r.h.s. of Eq. (gT]), 



T^f = Pf^f + ^--TiPf+Pf)Vf 



(S) 



Pi 



a , 
-f 3- I 
a 



-UJf)' 



k 



(29) 



(where ujf = Pf/pf is the equation of state parameter for the f-th fluid) corresponds to the gauge-invariant total 
(magnetic- field-induced) /-th fluid energy density perturbation [3a |. 
To derive Eq. (I7f|) we use the Einstein equation for 5tq, which is 



-Pf)vf^/k. 



(30) 



As a consequence of the high conductivity, cr 3> 1, of the primordial plasma, the Oi component of the magnetic field 
stress-energy tensor, cx E • B, being suppressed by I/cr, does not contribute to the r.h.s. of Eq. ^U\i . The energy 
density of the electric field, (x E'^, is suppressed by and does not contribute to the r.h.s. of Eq. ([?7|) . 

The trace of the space-space part of the Einstein equation gives an additional constraint equation. Since this 
equation is for the trace [131 , it does not have a contribution from magnetic anisotropic stress (but isotropic 

pressure does contribute), 



fc2 

* + — (* + $) 
o 



-AnGa^ I '^sjPf^f 




(31) 



Here c| ^ = dpf/dpf is the square of the speed of sound in the /-th fluid. 

Equations ^7} . ^E^ . and (pij) govern the evolution of the scalar metric perturbations $ and 4*, if the density, 
velocity, and anisotropic stress perturbations for each /-th component are known. For unmagnetized perturbations, 
the energy density and anisotropic stress of the magnetic field, pB(k, ry) and n('^^(k, 77), vanish. In this case Eq. ([^ 
results in the simple relation $ = The presence of collisionless particles, such as neutrinos, induces anisotropic 
stress [2^ . Even through we neglect neutrinos, in our case a stochastic magnetic field induces anisotropic stress 
and so violates the condition $ = — ^P. 

Using Eq. (^5]). setting pf, = = pc, and neglecting neutrinos, Eq. (PT|) becomes 



a V 



$ + - I 2$ - * 
a 



^ -2^ 
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Pi 
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5-y — 



3i^\ 
Pi ) 



(32) 



where the last step uses the zeroth order Friedniann equation in the radiation dominated model. Equations (|27p and 
(f28l) may be combined together, 



fc2(2$ + 5f) =47rGa^ 



Pi ( ^1 



3L(^\ 
Pi ) 



- PbSb + PcSc 



f 



(33) 



We show below that the combination 5^p^ — SL'"^^ on the r.h.s. of this equation also appears in the equations for 
matter perturbations and reflects an effective temperature rescaling (induced by the presence of the Lorentz force). 
Using Eq. HI]), Eqs. ([27]) and ([281) may be rewritten as 



f f 



(34) 



(35) 



These equations show that we should take note of two effects, one a consequence of anisotropic stress, the other arising 
from non-zero energy density and peculiar velocity perturbations. 



B. Matter perturbations 

The first-order energy conservation equations for photons, baryons, and the CDM fluid are 

^7 + IK^^ +4$ = 0, (36) 

Sb + fcwp' + 3$ = 0, (37) 
6c + fcw(^' + 3$ = 0. (38) 

The Lorentz force directly modifies only the Euler equation for the baryons, since only baryons are charged. The 
Euler equation for the CDM fluid is identical to that in a model without a magnetic field. Prior to decoupling, photons 
are tightly coupled to baryons and they move together. Since the Lorentz force affects the motion of baryons, the 
presence of a cosmological magnetic field also influences the evolution of photons. The first-order Euler or momentum 
conservation equations in the tight coupling regime (prior to last scattering) for photons, baryons, and CDM are, 

€^ - -k^ + f(«(^) - .p)) + f e(^) = 0, (39) 



-7 ^ ' ^"7 "b 



kL(s) 



K ' + -K ' - fc* - l5i<' - < ') + = 0, (40) 

a rt 



Pb 

yiS) + %(JS) _ fc^ ^ 0. (41) 

a 

In Eq. ((39|) Q^'' = 5n^y'^''/12 is the quadrupole moment of the photon temperature fluctuation and reflects the 
anisotropic nature of Thomson-Compton scattering. The photon density fluctuation is related to the temperature 
monopole moment, 6-y — AQq^\ and the perturbed photon velocity is the dipole term v^^"^ = q\^\ In Eqs. pop and 
()40p T is the differential visibility function; f = n^XeO'TO, where ne{z) is the charged particle number density, Xe{z) 
is the plasma ionization fraction, and ctt is the Thompson cross-section. In Eq. ()40p the term cx f (the so called 
"baryon drag force" term) reflects the coupling between photons and baryons, and so determines the velocity difference 

between the photon and baryon fluids. In the lowest order of the tight coupling approximation the t{v\ ) 

iS) _ JS) 



terms in Eqs. (|39p and (|3D]) vanish since 

At early time the hydrodynamical description for photons is a reasonable approximation due to their strong in- 
teraction with baryons and short mean free path. In this case II^'^-* = and so Oj"^' ~ 0. Subtracting Eq. ((40|) 
from Eq. pop . multiplying by R, and using Eq. (|37p. the velocity difference between the photon and baryon fluids, 

Au-yf, = ti^"^-* — u^'^'' , obeys 

L /Qr(S) \ DA 

RM^b + (1 + R)i/^v^b = T +R5^]- —{5b + 3$), (42) 

4 \ J ka 



where l^v-^h = w 



7 "b ■ 



i? cx a, SO at early times it is very small, and at these times Eq. 



results in"* 



Av^b 



(44) 



For modes with wavelengths larger than the Hubble radius, the difference between photon and baryon fluid velocities 
is significantly suppressed by the k/r factor in this equation. 



To derive an equation that describes the photon-baryon fluid in the tight coupling approximation, where 



iS) 



,{S) 



we multiply Eq. (gO]) by R and add Eq. ^ to get 



d_ 

drj 



]+{l + R)k^. 

Pi J 



(45) 



Here we have made use of the relation a/ a — R/R. Equation (j45p is valid over a limited range of times prior to 
decoupling. After decoupling photon evolution must be described by the Boltzmann transport equation and r is not 
large enough to ensure the equality of the photon and baryon fluid velocities, i.e., v^^^ ^ ''^b^\ so Eq. (|45p is not 
valid. Also, after decoupling it is no longer possible to ignore the photon anisotropic stress term (i.e., the temperature 



quadrupole moment H^'^'' = 1292''V5) that appears in Eq 



Using = 4ef ^ and = ep\ 



Eqs. ([36]) and (|45|) may be expressed in terms of temperature multipoles. 
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(46) 
(47) 



The quantity 0q + ^' is usually called the effective temperature, and in a model with a cosmological magnetic field 
it obeys the second order differential equation 
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drj 
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i?)(ef) 



iS) 
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(48) 



This equation is derived from Eqs. ([46]) and (l47|) . In addition to the usual baryon drag force term (x Rk^^/3 and 
gravitational potential time derivative difference term (x (vP — $), there is a term on the r.h.s. of this equation, 
oc k'^L^^^ / Pry, which directly reflects the presence of the cosmological magnetic field. 

Since L*^'^' iv)/ P-yiv) is time independent, we may use Oq'^' +'5 — 3L^'^-'/(4p7) as the generalized effective temperature 
in the case when a cosmological magnetic field is present. Equation 



~ ( S) 

rewritten in terms of 0X , where 40, 
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(49) 



Gq"^-* reflects the rescaling of the photon fluid energy density perturbation in the presence of a magnetic fleld. 



In a model without a cosmological magnetic fleld, deflning rrics 
(83)ofRef. [Mi), 



1 + R, Eq. 



can be rewritten as (also see Eq. 



3 ^" 



^0 



drj ( 



(50) 



For a model without a magnetic field Av-^i, obeys 

RAv^b + (1 + R)-i-Av^i, = R. 



kS^ a , ■ 

(5t, + 31.) 

4 fea 



(43) 



and so for R <g 1 Au^yj, vanishes. In the presence of a magnetic field the Lorentz force term is responsible for a non-zero velocity 
difference at early times. 



Equation (|50p is the second order differential equation that governs the dynamics of photon density perturbations 

= AQq^^). In the absence of gravitational potentials the r.h.s. of Eq. ((501 vanishes. The l.h.s. of Eq. ([SOI differs 
from the equation for an undriven simple harmonic oscillator only by the time-dependent factor rrieff. Just like the 
case for a harmonic oscillator, Eq. ([50|l has two independent solutions — sine and cosine modes — that depend on 
initial conditions. Defining the photon-baryon fluid sound speed cs — l/y/SmeS, for mcs/n^cS ^ w where uj — csk is 
the oscillation frequency, the JWKB solutions of Eqs. and ((17|) are, [3]| . 

e^^^ = AimJ^\osiks + <j)), (51) 
= yli%/3 m;ff^/''sin(fcs + 0). (52) 

Here Ai is the amplitude, s = J csdrj is the acoustic Hubble radius, and 4> is the phase. The constants Ai and (j) 
depend on initial conditions. 

We note that baryon pressure has been neglected, ph — 0, in the baryon Euler equation dJO]). Consequently Eq. 
also assumes that the baryon pressure vanishes. We discuss this assumption in the following subsection C. Here we 
assume vanishing baryon pressure and study acoustic oscillations driven by a weak Lorentz force. We again neglect 
gravitational potentials but now retain the last, Lorentz force, term on the r.h.s. of Eq. (|48p . In this case the leading 
JWKB terms in the equation are (again under the assumption that moff/'rioff ^ 

e<« + *^4ei"^|^^|4*^^. (53, 

The solutions of Eq. ([53]) are of a similar oscillatory form to those in Eqs. (|5ip and ((52|) . but now there is a constant 
shift of 8q'^^ Qq^^ + SpB/{4,py) (here we have used L^''^^ ~ —ps), while Q[^^ — — SOg'^^/A: remains unchanged. 

C. Acoustic oscillations in the baryon fluid 

The propagation of magnetosonic waves and magnetohydrodynamical instabilities in the expanding Universe are 
discussed in detail in Ref . [33| . The effects on CMB temperature anisotropics of magnetosonic waves in a homogeneous 
magnetic field are studied in Ref. [23|, where it is shown that a homogeneous magnetic field induces three types of 
MHD waves (fast and slow magnetosonic and Alfven waves) in an expanding Universe. 

Fast magnetosonic waves result in a rescaling of the fluid sound speed, i.e., cs —> \/c| + i;^ , where the original fluid 
sound speed cs is characterized by the fluid pressure p and energy density p, the Alfven speed va = Bq/ y^iir^p + p), 
and Bq is the (unperturbed background) homogeneous magnetic field strength ^] . Fast magnetosonic waves require 
a small inhomogeneous magnetic field Bi, so we write the total magnetic field B = Bq + Bi, where |Bi| <C |Bo|. The 
induction law in this case is 

^=Vx[vxBo], (54) 

where we work the leading order in the inhomogeneity. For the case of a homogeneous magnetic field, Adams et 
al. [l^] consider a zeroth-order (background) magnetic field Bq with zeroth-order energy density, pB oc Bq, small 
compared to the energy density of the photon-baryon fluid. In this case Bi is a first order perturbation. The fluid 
3- velocity perturbation v is also first order and satisfies the linearized Euler equation, [H, [sj , 

p|^+Vp+^[Box(VxBi)]-0. (55) 

Here the pressure gradient Vp is related to the sound speed through Vp — Cg'^p- In addition, the perturbed magnetic 
field obeys the Gauss law V • Bi = 0. 

Assuming tight coupling between photons and baryons, multiplying the baryon Euler equation by R and adding 
the photon Euler equation^ , Adams et al. [131 obtain the Euler equation for the photon-baryon fluid accounting for 
a zeroth-order spatially homogeneous background magnetic field (also see Eq. (52) of Ref. [39l|). 

^ [(1 + R)v] ~ cl,k6, = ^6., + (1 + R)k^ + ^-^-A__[k . {Bo X (k X Bi)}], (56) 



^ The baryon and photon Euler equations are Eqs. (11) and (13) of Ref. ]2^ . 



where cs.t is the baryon fluid (not the photon-baryon fluid) sound speed. 

Compared to the stochastic magnetic field case of Eq. , Eq. — for a homogeneous background magnetic 
field — contains an additional baryon pressure term on the l.h.s., Cg f^kSi,, while the term on the r.h.s. of Eq. (|56|) 

oc k • {Bq X (k X Bi)} is the analog of the Lorentz force term (x i'^"^-' on the r.h.s. of Eq. (|l5)) . We emphasize that 
Eq. (j45p is valid for a stochastic magnetic field and that here the smoothed amplitude of the magnetic field B\ is 1/2 
order in the perturbation expansion, while pb(cx i?^) and L*-'^' are first order. We argue below that at linear order the 
additional baryon pressure term (c| i,k5b) in Eq. ([55]) can be neglected. However, we emphasize that even if the sound 
speed in the uncoupled baryon fluid vanishes, i.e., cs^b = 0, the effective sound speed in the coupled photon-baryon 
fluid is not the same as the sound speed in the uncoupled photon fluid. Due to the tight coupling between photons 
and baryons the photon-baryon fluid sound speed depends on the baryon fraction, cs = l/v^Sm^ (see Eq. (|50p). so 
the coupling between baryons and photons reduces the sound speed from the l/\/3 value for an uncoupled photon 
fluid. 

The last term on the r.h.s. of Eq. ((5^ induces fast magnetosonic waves in the presence of a homogeneous magnetic 
fleld. These magnetosonic waves change the photon-baryon fluid sound speed, increasing it relative to the case without 
a magnetic field. In the limit of a weak magnetic field, the effective sound speed is cs = y/l/{3mes) + v\ [l^]. The 
Alfven speed va here is that defined in the photon-baryon fiuid, i.e., 

2 ^ Bl ^ ?,Bl 

^ 47r(pB + 4p^/3) 167r(l + i?)pT,' ^ ' 

For the more realistic case of a stochastic magnetic field, the Alfven speed should be defined in the terms of the 
smoothed magnetic field B\ [l^. Equation (|45|) is the analog of Eq. ([56]) for the case of a stochastic magnetic field. 
For this stochastic field case we define the Alfven speed squared as v\ = 3pb/{2{1 + R)pj). Here we neglect the 
baryon pressure, pf, ~ 0. Note that for the case ns — —3 our definition of the Alfven speed coincides with that 
used in Ref. [2J| under the assumption that Bx = Bq, i.e., va = va- The term cx L^"^^ on the r.h.s. of Eq. 
ensures that the effective sound speed is rescaled in a manner similar to that for an homogeneous magnetic field, 
Cs -^cs = y/l/{3mcs) + v\. 

This rescaling of the sound speed may be formally described as a baryon energy density fraction change — > i? = 
R — AR, as follows. We define fhcs — c|/3, then 

™eff 

"^'=« = TT9Bf7(I^' ^''^ 

so the rescaling of the sound speed induced by the presence of a magnetic field (a homogeneous or stochastic field) 
is equivalent to a reduction of the baryon fraction, AR = 3w^m^jj/(l -I- 3v\mcs)- Consequently, this increase of the 
sound speed (relative to that of a model without a magnetic field) induces shifts of the CMB anisotropy angular power 
spectrum peaks comparable to shifts resulting from a reduction of the baryon density pb —> pb — 3{l + R)BQ/{ATr); here 
we have used the fact that va cs- This was noted, for a homogeneous magnetic field, from the result of numerical 
simulations, in Ref. [131; the analytical results we have derived, for a homogeneous or a stochastic magnetic field, are 
new. The reduction of the baryon fraction reduces the baryon drag force (cx Rk^'if) on the r.h.s. of Eq. As a 

result there are two effects: a shift of the CMB anisotropy angular power spectrum peak positions and a reduction of 
the peak amplitudes. 

We also note that the Lorentz force in the Poisson equation ([34]) can be treated as reducing the photon effective 
temperature. In particular, photon energy density perturbations {S^) can be compensated by the Lorentz force if 
6^{ri) = 3L^^^ / p-y (also see Eq. ((55|) ). Note that L^^'> / p-^ is time independent. 

Ref. [IB'l (their Eq. (46)) considers a modified form of Eq. (j56p for the case of a stochastic magnetic field: they 
discard the Lorentz force term. Ref. ^] also neglects the Lorentz force contribution (their Eq. (13) and App. A). 
Both references argue that such a force-free approximation is justified by the infinite conductivity of the plasma which 
results in a vanishing electric field in the metric perturbation equation. We have noted however, at the end of Sec. 
II, that the force-free approximation cannot be used in this manner. Here we point out that the current, cx V x B, 
does not vanish, and so a Lorentz force term must be present in the baryon Euler equation. 

We have shown that the Lorentz force term in the baryon Euler equation results in shifts of the CMB anisotropy 
angular power spectrum peaks, relative to the case without a magnetic field. ^ This result, obtained for a stochastic 
magnetic field, generalizes that of Ref. [13] for an homogeneous magnetic field; for = —3 our estimate reproduces 



Below we discuss the peak shifts related to the baryonic pressure effect, which are small compared the shifts discussed here. 



the acoustic oscillations shown in Fig. 1 of Ref. [24]. On the other hand this result contradicts Fig. 1 of Ref. [27j . 
This figure is for the case of a stochastic magnetic field and it indicates that the magnetic field effectively increases the 
baryon fraction (while Ref. |27i] considers a magnetic field with positive spectral index, ng > 0, this cannot explain 
the effective increase of the baryon fraction they find) . 

We note that Ref. retains the baryon pressure term c>c c| ^kSb (and argues that c| ^ is related to the gradient of 
the magnetic pressure, Vps) in the photon-baryon Euler equation Ref. !26] claims that their CMB anisotropy 

angular power spectrum peak shifts are similar to those found in Ref. 24] , because of this term. The Cg i^kSb term 
does induce a rescaling of the sound speed. However, this term is second order {St and Vps are first order), does not 
contribute at linear order, and so can be discarded compared to the Lorentz force contribution which they neglect. 



D. Initial conditions 

(S) 

The presence of a cosmological magnetic field modifies the initial conditions for the monopole 8q and gravitational 
potential, $ and 4", perturbations. A proper treatment in an inflation model requires analysis of quantum mechanical 
fluctuations during inflation, see, e.g., Refs. [13] for the case without a magnetic field. Here we adopt a more 
phenomenological approach. 

As discussed in Refs. 0, [13], there are three different types of perturbations to consider, depending on conditions 
in the radiation dominated epoch: (i) adiabatic, where the initial gravitational potentials are large compared to the 
magnetic field fraction of the energy density, i.e., ps/p-y <C ^in and ps/p-i <C $in, so the magnetic field energy density 
may be ignored; (ii) quasi-adiabatic, where ps/p-y < ^in, Ps/p-y < ^in', and, (iii) isocurvature, where the magnetic 
field energy density fraction dominates over the initial gravitational potentials, ps/p-y 3> ^in and ps/Py 3> $in- 

Since adiabatic and quasi-adiabatic CMB perturbations have been studied in some detail (for a review see Ref. 
[39|), we focus on isocurvature fluctuations induced by a cosmological magnetic field. 

We assume (as is conventional in the case of an isocurvature solution) that initial values of all relevant variables are 
determined by the magnetic fleld energy density and anisotropic stress. Under such an assumption the initial conditions 
for the gravitational potentials are obtained through Eqs. (|27p and (|35p (assuming that initial fluid perturbations are 
zero), 

,2^ 47rG , , 
k <I>in = -^PB[m), 

k\^,,, + ^,,)^-^UiS){rjo), (59) 
< 

where Oin is the value of the scale factor when the initial conditions are applied. Adding Eqs. ([5^ and using Eq. (fTil) 
implies that fc2(2$i„ + ^-in) = -12T:GL'^'^\r]o)/af^. 

Using these initial conditions for the gravitational potentials and assuming a weak magnetic field, one may obtain 
solutions for scalar magnetic perturbations with wavelengths larger than the Hubble radius (i.e., the leading terms in 
an expansion in ki] <^ 1). 



IV. CMB TEMPERATURE ANISOTROPIES 



In this section we compute the CMB temperature anisotropics due to the presence of a cosmological magnetic field. 
We assume the existence of a cosmological magnetic field on scales larger than Hubble radius, thus we assume that 
a magnetic field has been generated during inflation Q. Our analysis below holds for a magnetic fleld with spectral 
index ub larger than —3. 

CMB temperature fluctuations are caused by scalar perturbations due to: i) initial intrinsic inhomogeneities on the 
last scattering surface; ii) the relativistic Doppler effect due to the baryon velocity as the photon propagates to the 
observer; iii) the difference in the gravitational potential between the points of photon emission and the observer (the 
usual Sachs- Wolfe effect, SW); and, iv) changes in the gravitational potential as the photon propagates (the so-called 
integrated Sachs- Wolfe effect, ISW) P, [H, |4l| . 

Using the total angular momentum formalism (3lj ] , the angular power spectrum of the CMB temperature anisotropy 
measured today is 

-^r 2£T^- 



Here I is the multipole index and (vo, k) is the l-th multipole moment of the (scalar-sourced) temperature fluctu- 
ation and is determined by the integral solution of the Boltzmann transport equation [3l| , 

(61) 

Here P^'^^ = (62'^'' — \/6-£'2'^'')/10 is the anisotropic (quadrupolar) part of the Compton scattering cross-section 
e of CMB p 

quadrupole moments^, and the radial functions 



(5) (5) 

which is a source of CMB polarization anisotropics — where Q2 ^^'^ ^2 temperature and ^^-polarization 



= 3d^), = j'Ax). 3f'\x) - \wi+n[x)l (62) 

where ji is the spherical Bessel function and a prime represents a derivative with respect to x. 

Equation ([61]l includes the four effects responsible for the CMB temperature anisotropics mentioned above. The 
initial photon temperature (oc P^^'^), SW (oc (Og'^-* -I- \1/)), and baryon velocity Doppler (oc w^'') effects are present on 
the last-scattering surface and so in Eq. (|6ip they appear with the factor e^'^f; the ISW effect, cx (4' — $), contributes 
from decoupling until today and thus is suppressed by the factor e""^. 

We are interested in the large-scale CMB temperature anisotropy due to a cosmological magnetic field. The 
contribution from the ISW effect is negligible when compared to the SW effect. Also, compared to the SW effect, 
the Doppler term oc v^^^^ plays a secondary role when 1 + i? > 1 [3l|. The quadrupole term P^'^^ oc kQ^f^ /f (see Eq. 
(90) of Ref. ^ and the discussion around Eq. (I75|) in Sec. V below) is strongly suppressed for kr] <^1, and thus we 
neglect it on large angular scales. On large angular scales the largest contribution arises from the ordinary SW effect 
so we approximate the temperature integral solution for the scalar perturbation as 

^2^^ ^ l\ e-^f(e|,^) + ^)Mkvo - M- (63) 



The visibility function fe is sharply peaked at decoupling, so we use the approximation, 41 1 



''21+1 - [Qo''^(^doc) + ^{Vdcc)]je{kvo - kmcc) ^ -^^mccMkvo), (64) 

where Tydcc is the value of conformal time at decoupling. This is the familiar SW, $/3, result. This expression for the 
temperature anisotropy multipole moment depends on <i>(?7dcc), so we must solve for $(?]). 
The gravitational potential <l>(r/) obeys Eq. ((?71) . To solve this equation we decompose $ as 

$(77) = $i(77) + $2(r?), (65) 

where the potential $1 is due to the magnetic field energy density and <I>2 is related to the energy density and velocity 
perturbations in the fluids (these perturbations in the fluid are induced by the magnetic field anisotropic stress) . From 
Eq. ([27]), the potentials ^i{r]) and $2(??) obey 

fc2$i = AnGa^pB, (66) 
/c2$2 = AnGa^^Vf, (67) 
/ 

where 2?/ is the gauge- invariant energy density perturbation in the f-th fiuid, Eq. (|29p . 

Mathematically Vj should be obtained through the solutions for the energy density {5f) and velocity {v^ ) pertur- 

bations of the f-th fluid. Sf and obey fairly complicated second order differential equations (see Sec. HI.B) that 
are not straightforward to integrate. However, since we only consider perturbations arising from a magnetic field. 



The zeroth-order term in the expansion of the photon Boltzmann transport equation for £ = 2 relates ©1^' and si^^ as si^^ 



-V60'^V4, leading to Pf^) = g!] 



2 ""^ -^2 

^ ) /A mil 



it is expected that on scales larger than the Hubble radius the sum of the density and velocity perturbations in the 
fluids, 2?/, should be of order the magnetic field energy density, while on small scales radiation pressure prevents 
perturbations from growing [l^ . Thus X)/-^/ — Pb-^ Using this, Eqs. (|66p and (|67p imply that the maximal value 
for $ is 

With this expression for the temperature anisotropy multipole moment, Eq. (j64p . becomes 

ef^(ryo,fc) 8^G 



PBivo,k)ji{krio), (69) 



where adoc is the value of the scale factor at decoupling. Using this in Eq. ((5n|) . the CMB temperature anisotropy 
angular power spectrum is given by 

cr'-"^§(|g)7°°..^%ii^,i»,.)^ (70, 



Now \pb\^ is given in Eq. so in terms of the Bessel function — y/2x/Tr ji{x), we find 

JhiMkVo). (71) 
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The integral in Eq. ([7T|) may be evaluated by using Eq. (6.574.2) of Ref. [42]. For ns < —3/2, when the magnetic 
source is dominated by the term proportional to fc^"'^^^, we find 



«2^ee(S) _ 22"«+inBr(-2ng) (5) 2„b+2 ^7o^ 
-(nB+3)r2(-nB + l/2)^««^ ' ^^^^ 



while for ub > —3/2, when the magnetic source is dominated by fc^"^"'"'^, we have 



In these expressions 

^(5) ^ ^2g2j^4^2nB+6 

®® 3(2ns + 3)r2(nB/2 + 3/2)a4^,r;2"«+2- 



Equations (j72p and (j73p describe the angular power spectrum of the CMB temperature anisotropy induced by scalar 
magnetic perturbations. The maximum growth rate of the power spectrum ("^Ce with £ is i^-^. This occurs for 
Ub > —3/2. So, as expected, the CMB scalar temperature fluctuations due to a stochastic magnetic field are strongly 
suppressed at small angular scales. The suppression is stronger for an inflation generated magnetic field with ub ^ —3 



V. CMB POLARIZATION ANISOTROPY 



In this section we compute the scalar magnetic-field-induced CMB i?-polarization anisotropy angular power spec- 
trum. Scalar (density) perturbations only induce electric type i?-polarization anisotropics. 



We justify this approximation by using the fact that the evolution of the magnetic-field- induced energy density and velocity perturbations 
on large scales must follow the dynamics of the magnetic source. On the other hand, the homogeneous solutions of Eq. H27|l on small 
scales (in a model without a magnetic field) are sound waves (oscillating energy density perturbations). To obtain solutions of the 
inhomogeneous equation we must integrate the product of the source and the homogeneous solutions; as a consequence of the oscillation, 
the resulting integral will be negligibly small. 



In the total angular momentum formalism 31], the £^-polarization anisotropy angular power spectrum is 



where Ei is the l-th multipole moment of the i?-polarization anisotropy. It is determined by the integral solution of 
the Boltzmann transport equation [sij, 

rvo 

^ = -\/6 / dT^ e-^T F(^)(77) ef ^(fcr;o - fcry). (76) 



2^+1 

Here P^^^ is defined below Eq. (pT|) and e^'^'' are the i?-polarization radial functions 



Polarization anisotropics, being connected with the anisotropic stress, are determined by shear viscosity. They are 
generated during last scattering when the tight coupling is enhanced by the fast growth of the inverse differential 
visibility function l/f [3l|. In the lowest order of the k/f expansion, the set of Boltzmann equation solutions are 
given in Eq. (90) of Ref. [31[, and for the scalar perturbation mode, 

9r 

Here 8^'^'' is the dipole moment of the temperature anisotropy and obeys Eq. ((17|) . A similar equation holds for the 
case of scalar magnetic perturbations, see the discussion in Sec. 3 of the last of Refs. ^19]. In this case, in addition 
to fluid shear viscosity, there is also a contribution from magnetic anisotropic stress Since we focus only on 

magnetic-field-induced effects, we consider H^'^^ to be the dominant source of polarization. Then, according to the 
technique used in Sees. IV. B — D of Ref. [3l|, and using footnotes 2 and 7, we have 

P^'Hv) = ~-^nis^-^PBivo,k). (79) 
Using these expressions for P^^'' and e^'^', Eqs. (f7^ and d77|) . in Eq. ([TB]) . we find 



0^,1- o /7f / 1 ^jfXkiio- krj). (80) 

Here we use again the fact that the visibility function e^'^r peaks at decoupling and approximate the ii'-polarization 
integral solution as 

^ -Je(kr]o). (»1) 



2^+1 ■iV2p'yoT{T]dcc)Vo 

To obtain the angular power spectrum cf^^^^ we need T(?7doc)- From Eq. (C3) of Ref. [4l|, and assuming for the 
current value of the baryon energy density parameter rif,(?7o) = pb/ Per = 0.05 and that the redshift at decoupling 
Zdec = 1100, we get 

f(zdec) = 8.05 (-\ . (82) 

V«/doc 

Photon-baryon decoupling occurs during the matter dominated epoch, when a(ri) cx 77^, so (a/a)doc — '^/Vdcc, and 

f{r]dec) ^ 16.1/77d 

Using Eqs. ([HT|) . (IH^ . and the CMB _E-polarization angular power spectrum of Eq. (fTS]) is 



3 X 2i4(2nB + 3)r2(nB/2 + 3/2)^2^ ^ ^ 



^ n_B + 3 



JhiMkVo). (83) 



To evaluate this integral we consider two cases: ns < —3/2, when the integral is dominated by the second term 
(oc and, ub > —3/2, when the main contribution to the integral comes from the term cx fc^"^'^'^. In the first 

case, for —3 < Ub < —2, the integral may be exactly evaluated using Eq. (6.574.2) of Ref. [i^).^ 



e2(jEE{S) UBThnB - 2) (S) 2ns + W (o < < _2) f84) 

- Av^iuB+m-nB-m '''' { d<nB< 2). m 



For —2 < hb < —3/2 we may evaluate integral using the semi-analytical approximation of the Appendix of Ref. [l8|. 
For X 3> 1 Ji+i/2{^) — ■\/2/ (ttx) cos[x ~ {I + l)7r/2], Eq. (9.2.1) of Ref. [i^, and replacing the oscillatory function 
cos^ X by its r.m.s. value of 1/2, we get (see Eq. (B2) of Ref. [3]) for ub = —2, 

^2c^msj ^ _2 In Afie^ [UB = -2), (85) 

while for -2 < ns < -3/2, 

^'^^""^'^ =^ 9( ^o^ (fc^^o)^""+^4''i^^ (-2 <nB< -3/2). (86) 

Z(nB + 'ijyriB + Z) 

For the case when ub > —3/2, using x^j^'^^^ J^^" dx x Jf^i/2i^) — ^^u/tt, we have 

^2cEEiS) ^ 8 ^^^^^^2„.+4^(S)^6 (^^ > _3/2). (87) 



In these expressions 



.(5) BtX^-^+^rjl 



' — 



3 X 2i4(2nB + 3)r2(ns/2 + 3/2)fjl''-+'° p^^^ " 

Contrary to the temperature anisotropy, the i?-polarization anisotropy angular power spectrum i'^Cf^ grows rapidly 
with £ (the fastest growth rate is £^). An £^ dependence for i?-polarization is also expected in the absence of a 
magnetic field [3l| . 



VI. TEMPERATURE— POLARIZATION CROSS CORRELATIONS 



In this section we compute the scalar magnetic-field-induced CMB temperature — £'-polarization cross-correlation 
anisotropy angular power spectrum [3l| . 

- ^ / ^ 2£ + i 2eT)~- ^ ' 



Using Eqs. (HH), dMl), fHT I) . and ((5^ . we find 

Ji+i/2{km), (90) 



3 X 210(2^5 + 3)r2(nB/2 + 3/2)7^^al^^p^o J k 



^ 775+3 



under the approximations discussed in the two previous sections. 
Evaluating the integral in Eq. we find, for —3 < ub < —3/2, 



^{nB + 3)T{-nB- 1/2)' 



while for ub > —3/2, 



The magnetic field source is non-zero up to the damping scale k u . Since the integral is dominated by small wavenumbers we can replace 
the upper-cut-off scale k]j by oo. 



where 



3x2ii(2nB + 3)r2(nB/2 + 3/2)a2^^ryg"«+V^o- ^ ^ 



As in the case of temperature fluctuations, the temperature — i5-polarization cross-correlation angular power spectrum 
I'^Cf^'^^^ has maximum growth rate c>c for ub > —3/2. 



VII. CONCLUSIONS 



We present a systematic discussion of scalar isocurvature magnetic perturbations (magnetosonic cosmological waves) 
in a Universe with a stochastic primordial magnetic field. We derive the complete set of equations that govern the 
dynamics of linear magnetic energy density perturbations. 

A stochastic magnetic field shifts the acoustic peaks of the CMB temperature anisotropy angular power spectrum, 
acting in a similar way as a reduction of the baryon fraction. This result extends the work of Adams et al. [23| who 
studied a homogeneous magnetic field. The second important effect that comes from a stochastic magnetic field is a 
non-zero anisotropic stress which generates a CMB i?-polarization anisotropy. 

We obtain approximate analytical expressions for the CMB anisotropy angular power spectra cf^^^\ Cf^'^^\ and 
^QE(s) ^ Numerical values of these spectra depend on four parameters: the cut-off wavenumber ko', the smoothing 
length A; the amplitude of the smoothed magnetic field By, and, the magnetic field power spectral index ub- 

We find that the scalar CMB temperature anisotropy power spectrum amplitude cf^^^'' rapidly decreases with 
increasing consequently, the magnetic field energy density contribution to the total CMB temperature anisotropy 
signal is suppressed at large multipole number. On the other hand, the contribution from the magnetic field anisotropic 
stress to the i?-polarization anisotropy becomes large on small angular scales and should be accounted for when 
estimating the magnetic field contribution to the CMB £'-polarization anisotropy (it must be added to the vector 
CMB £;-polarization) [13, [ll]. 

Since scalar CMB temperature perturbations induced by a stochastic magnetic field do not depend on magnetic 
helicity, precise measurements of the CMB temperature anisotropy angular power spectra peak positions and am- 
plitudes, combined with a CMB-independent measurement of the baryon fraction and CMB polarization Faraday 
rotation data, should provide information on the symmetric part of the magnetic field power spectrum, Pb- This 
determination of Pb together with future CMB _B-polarization anisotropy data (which is sensitive to both power 
spectra, Pb and Ph), could lead to a constraint on cosmological magnetic helicity. 
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